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Abstract 

It is presented a theory that describes a spin glass phase at fi- 
nite temperatures in Kondo lattice systems with an additional RKKY 
interaction represented by long range, random couplings among local- 
ized spins like in the Sherrington- Kirkpatrick (SK) spin glass model. 
The problem is studied within the functional integral formalism where 
the spin operators are represented by bilinear combinations of fermionic 
(anticommuting) Grassmann variables. The Kondo and spin glass tran- 
sitions are both described with the mean field like static ansatz that 
reproduces good results in the two well known limits. At high temper- 
atures and low values of the Kondo coupling there is a paramagnetic 
(disordered) phase with vanishing Kondo and spin glass order param- 
eters. By lowering the temperature a second order transition line is 
found at T$g to a spin glass phase. For larger values of the Kondo 
coupling there is a second order transition line at roughly T^ to a 
Kondo ordered state. For T < T$g the transition between the Kondo 
and spin glass phases becomes first order. 



1. Introduction 

The antiferromagnetic s-f exchange coupling of conduction electrons to local- 
ized spins in heavy fermion rare-earth systems is responsible for two compet- 
ing effects: the screening of the localized moments due to the Kondo effect 
and the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction among mag- 
netic impurities which may induce a long-range magnetic (antiferromagnetic 
or ferromagnetic) ordering or eventually a spin glass magnetic ordering. The 
Doniach diagram[|lj gives a good description of this competition: the Neel 
temperature T/vis firstly increasing with increasing the absolute value of the 
exchange interaction constant Jk, then it is passing through a maximum and 
finally it tends to zero at the "quantum critical point" (QCP). Such a decrease 
of T/v down to the QCP has been observed in many Cerium compounds, such 
a.sCeAl 2 , CeAg or CeRti2Si 2 , under pressure. Above the QCP, there exists a 
very strong heavy fermion character, but several possible behaviours, i.e. the 
classical Fermi liquid one with eventually a reduced Kondo temperature |2|. [| 
or different Non-Fermi-Liquid (NFL) ones, have been observed in Cerium or 
Ytterbium compounds]!, ||. 

But, in the case of disordered Cerium alloys, the disorder can yield a 
Spin Glass (SG) phase in addition to the NFL behaviour at low tempera- 
tures around the QCP. The magnetic phase diagram of CeNi\- x Cu x has been 
extensively studied^, [F|. CeCu is antiferromagnetic below 3.5 K and CeNi 
is a non magnetic compound with an intermediate valence. The low temper- 
ature antiferromagnetic phase changes, around x = 0.8, to a ferromagnetic 
one which finally disappears around x = 0.2. At higher temperatures , a 



Spin Glass state is deduced from all measured bulk properties, such as for 
example the ac susceptibility; for example, for x = 0.6, the SG state exists 
between 2 K and the Curie temperature T c = 1.1K. At x = 0.2, there ex- 
ists below 6 K only a SG state which transforms to the intermediate valence 
CeNi, after passing probably through a Kondo lattice system. Thus, in the 
CeNii_ x Cu x system, there appears a SG phase at finite temperatures above 
the ferromagnetic order and finally a SG-Kondo transition occurs for small 
x values. 

Another disordered Cerium system, namely CeCoGe 3 _ x Si x alloys, has 
been also studied for different x values, by different experimental techniques 
including muon spin relaxation || ||. The compound CeCoGe^ is antiferro- 
magnetic below 21 K and CeCoSi^ is an intermediate valence compound ; 
the QCP of these alloys lies around x = 1.5. The muon spin relaxation exper- 
iments in the quantum critical region (x — 1.1 to 1.5) show that a fraction of 
Ce ions experience random f-f indirect exchange interactions, which causes 
frustration of some Ce spins in the system. So, for x = 1.1 and x = 1.2, 
frustrated moments of Cerium freeze like in a spin glass while the other Ce 
moments form a disordered antiferromagnetic state. Thus, near the QCP, a 
spin glass state can exist in these alloys, in addition to the observed NFL 
behaviour. So, a striking novel behaviour, i.e. the appearance of a spin glass 
state at finite temperatures in some disordered Cerium alloys, has been ob- 
served and a SG-Kondo transition is developing with increasing x around the 
QCP. 

The aim of our paper is to present a theoretical model that describes 
the spin glass-Kondo phase transitions, and that we do by studying a sys- 



tern Hamiltonian that couples the localized spins of a Kondo lattice with 
an additional long range random interaction, like in the Sherrington- Kirk- 
patrick spin glass model [HJ. A similar Hamiltonian has been considered 



in [J. 1] to describe NFL behaviour and a QCP in some heavy fermions com- 
pounds, although the relevant approximations differ in this work and ours. 
The authors in [0]] are primarily interested in the description of the QCP at 
T = 0, then they solve first for the Kondo effect by decoupling the conduc- 
tion electrons bath into independent conduction electron "reservoirs", with 
no communication between the reservoirs at different sites. 

During the course of our work, another paper using the same Hamiltonian 
has been proposed^]. The representation of Popov and Fedotov eliminates 
the unwanted spin states but the approximations involved are in fact different 
from those used in our work. They study essentially the spin glass state and 
they finally obtain a second order SG transition with a transition temperature 
depressed by the Kondo effect, in second order perturbation theory. 

In the present paper we take a different approach: the localized spins of 
the Kondo lattice will be effectively immersed in a common bath of conduc- 
tion electrons and the Kondo effect will be studied in a quantum static ap- 
proximation that is basically equivalent to the mean field decoupling scheme 



[y, [L3| . The spin operators are represented by bilinear combinations of 
fermionic creation and destruction operators, for the localized f-electrons, 
and the spin glass transition will be studied within the static approxima- 
tion. This deserves some special discussion. In the Ising quantum spin glass 



(QSG) model |14jl , the spin operator S* commutes with the particle number 



operator rii S and the static Ansatz gives the exact answer, as the problem is 



essentially classic. When we add to the fermionic Ising QSG a s-d exchange 
coupling of the localized spins to the conduction band electrons the problem 
ceases to be exactly soluble and the static Ansatz is just an approximation, 
that we consider justified to describe a transition at finite temperature. It 



has been shown in [15] that the exact numerical solution of Bray and Moore's 



equations |T6| gives for the spin-spin correlation function Q{r) roughly its 



constant classical value at finite temperature, what justifies the use of the 



static Ansatz of [[16] at not very low temperatures in the Heisenberg spin 
glass. 

We consider then the static Ansatz which corresponds to an approxima- 
tion similar to mean field theory, where by neglecting time fluctuations we 
can provide a description of the phase transitions occurring at finite temper- 
ature. 

We use functional integral techniques where the spin operators are rep- 
resented by bilinear combinations of fermionic (anticommuting) Grassmann 
fields. As we show in the next section, this method is ideally suited to de- 
scribe a Kondo lattice transition, and it has been recently applied by two of 



us to the study of fermionic Ising spin glasses with local BCS pairing [[17 
Recent work |18] also showed the existence of several characteristic tempera- 
tures in the Ising fermionic model, with the de Almeida-Thouless instability 



HJ occurring at a temperature lower than the spin glass transition temper- 
ature. 

This paper is organized as follows: in Sect. 2 we describe the model and 
relevant results, we reserve Sect. 3 for discussions and conclusion, while the 
detailed mathematical calculations are left for the Appendix. 



2. The model and results 

We consider a Kondo lattice system with localized spins Si at sites % — 1 • • • N, 
coupled to the electrons of the conduction band via a s-d exchange interac- 
tion. It is necessary to introduce explicitly the resultant RKKY interaction 
by means of a random, long range coupling among localized spins like in 
the Sherrington Kirkpatrick (SK) model for a spin glass. To describe the 
Kondo effect in a mean-field-like theory it is sufficient to keep only the spin- 
flip terms || in the exchange Hamiltonian, while the spin glass interaction 
is represented by the quantum Ising Hamiltonian where only interact the 
z-components of the localized spins |TT], [TJ], [T7J . 
The Hamiltonian of the model is: 

H - fi c N c - fi f N f = H k - fi c N c - fi f Nf + H SG (1) 

H k - fi c N c - fi f N f = ^ ^n k a + e J2 n L + JkJ2 l S fi s ci + S ji s ti] ( 2 ) 

k,ff i,o i 

KSG = ~ Yl J iJ S fi S fj ( 3 ) 



i -J 



where J k > 0, 



Sji = filfo ! s ci = 4i d H ( 4 ) 



1 



t f *t 



Sfi — qUiT^T — fiifiil 



2 
5 



and /j£, fi a (df a , di a ) are creation and destruction operators for electrons with 
spin projection a =| or \ in the localized (conduction) band, that satisfy the 
standard fermion anticommutation rules. We also have n^ = d\ a d^ a where: 

VN V 



C^Ee-*** (5) 



The energies eo(efc) are referred to the chemical potentials /i/(/i c ), respec- 
tively. 

The coupling J^ in eq.(|3|) is an independent random variable with the 
gaussian distribution 



PW ,)= e -W^_ (6) 

Functional integration techniques have proved to be a suitable approach 
to describe phase transitions in disordered quantum mechanical many-particle 
systems ||20|| . The static approximation within this formulation consists in 
neglecting time fluctuations of the order parameter, and when it is combined 
with the neglect of space fluctuations it leads to the usual Hartree-Fock, 
mean field like approximation. When dealing with the Hamiltonian in Eq. 
(|l])-Eq.(|3|), we notice that in the limiting case 4 = we obtain a pure quan- 
tum Ising spin glass where the static approximation gives the exact result 



14, 18], while for J = we recover the mean field approximation that has 



been used succesfully to describe the Kondo lattice || |13| . Then we consider 



that the use of the static approximation has an interpolation character and 
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will provide reliable results to describe critical behaviour at finite tempera- 
ture |15| in systems that do not present a quantum critical point. 



In the Lagrangian formulation [|17], |2(| the partition function is expressed 
as 

Z= J D^VX^)^ (7) 

where the action A is given by 

^ = E jf *r{( j^lto* + (^Vi)lM - jf «(r)e*r (8) 

In both expressions ipi a (T) and ^(t) are anticommuting, complex Grass- 
mann variables associated to the conduction and localized electrons fields, 
respectively, while r is an imaginary time and j3 the inverse absolute tem- 
perature. 

We show in the Appendix that in the static, mean field like approximation 
the action A may be written: 

A = A + A K + A SG (9) 

with 

^o = J2J2i( iuJ - (3eQ)8ijijj\ a {ixi)ijj ia {uj) + (iuSij - Ptij)<ptr(w)<Pj*(w)] (10) 
where from Eq.(|36|) 



A 



K 



N 



EEC( w W")]E?UVw(oO] 



(11) 



A-SG — E JijSi^j 



i -J 



and in the static approximation [14, 17 



(12) 



The sums are over fermion Matsubara frequencies u = (2n + l)n. 
The Kondo order is described by the complex order parameter 



(13) 



1 

TV 



AU-^E^tM^H) 



A ^ = ^E(vLM^H) 



(14) 



that in a mean field theory [|, [L3j describes the correlations Aj. = (fi a di a ) 
and A CT = (dj<r/i<r) - Complex conjugation of Grassmann variables is defined 
through the transposition rule ]20| (ip'tp)* = (p'lf). 

We show in the Appendix that standard manipulations give for the aver- 
aged free energy within a replica symmetric theory: 



(3F = 2(3 J k X 2 + -p 2 J 2 (t + 2qx) ~ PSl 



(15) 



where 

N 



I /-+oo N n r+oo 

pn = ]im—{ U D ^U n^; ex p(EM^(u;)|)-i} ( 16 ) 

and the order parameters q, x, and A must be taken at their saddle point 



value. Here q is the SG order parameter |I3, [T7| and the static susceptibility 
is x = /3%. We use the notation Dx = -j=e~z x . 

The function Gij aa (u) in Eq.(^) is the time Fourier transform of the 
Green's function Gij a (r) = i(T/j o -(r)/j (T (0)) for the localized electrons in the 
presence of random fields Zj and e a j at every site, and from Eq.(^) satisfies 
the equation 

Gr.i» = [iu - [3e - ah jc \b io - (3 2 J 2 k \ 2 1%3 (u) (17) 

where 

h ja = ^(3J Zj + ^2xf3Je aj (18) 

while jij(u)) is the time Fourier transform of the conduction electron Green's 
function 7y(r) = «(Tdj o .(r)d| o .(0)) and is given by 

% 3 l = [iu - PudSij - PUj (19) 

We obtained in Eq. ([jj]) the Green's function for the f-electrons in a 



Kondo lattice ||13|| , but now in the presence of a random field hj at every site 
that prevents us from proceeding with the calculation. In the pure SG limit 
Jk = the Green's function in Eq. (|T7|) is local and the integrals in Eq. (|16|) 
reduce to a one site problem, while in the Kondo limit J = the random 
fields vanish and the integrals separate in reciprocal space. We adopt here a 



9 



decoupling approximation that is reminiscent of the model with independent 



'reservoirs" of electrons considered in [pT |. We replace the Green's function 



Gijaaiyj, {hi a ...hj a ...h Na }) by the Green's functions r^ U(T (uj,hj a ),j = 1...N, 
of N independent Kondo lattices, each one with a "uniform" field hj a at 
every site /x, v, by means of the approximation: 

\n\G£ a Mh 1 ...h N })\ « ^J2^K»Mh ja )\ (20) 

3 

where Y ^^iu.hj^) is the f-electron Green's function for a fictitous Kondo 
lattice that has a uniform field hj a at every site /x, z/ and satisfies the equation: 

r ^<r( w > M = [«*> - Z^eo - (x/ija] ^ - (3 2 Jl\ 2 ^^ (21) 

where, from Eq. (|l9l): 

VM^E^- (22) 

Now Eq. ( pTD may be easily solved by a Fourier transformation with the 



result: 



ln|r;J>,M| =£ln[I7>,/^)] (23) 



where 



rrV, M = [«" - /&„ - <rh ja ] - (3 2 J 2 k \ 2 . l — . (24) 

iu — pe k 

We may now introduce Eq. ( ^0|) and Eq. (|23|) in Eq. (0), the integrals 

over the fields separate and we obtain 

/+oo r+oo | _, 

Dz\n{ Deex V (Y,TfT, S ^ h ))} ( 25 ) 

-oo J— oo — iV _. 
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with 

^(£,/0 = £Mr fe ->,/0] (26) 

and h is given in Eq. (|18|), with z and e in place of Zj and Sj a . 

The sum over the fermion frequencies is performed in the standard way 
by integrating in the complex plane p0| , with the result: 



S a (k, h) = In [(1 + e-^+)(l + e-^-)l (27) 

where 

u a ± = ^{(3e k + ah] ± { l -{f3e k - ah) 2 + ((3J k \) 2 }l (28) 

We consider eo = that corresponds to an average occupation (tif) = 1, per 
site. 

Replacing sums by integrals , in the approximation of a constant density 
of states for the conduction band electrons, p(e) = p = -^ for —D < e < D, 
we obtain from Eq. ( fffi) to Eq. (J27D the final expression for the free energy 
in Eq. flRj|): 

1 r+co r+oc 

(3F = 2(3 J k X 2 + -f3 2 J 2 (x 2 + 2<?x) - / Dzln { / 7Jee E(/l) } (29) 
with 

£(/>) = — / cfclnlcosh l ; + cosh(\/A)}, (30) 

P-D 7-/3D 2 

A = ^(a;-/i) 2 + (/3J fe A) 2 (31) 

and from Eq. (|TJ|) we have h = f3J[^/2qz+\/2%e]. The saddle point equations 
for the SG order parameters are: 

" = L d ^-JdJ" } (32) 

11 



* = L Dz Td^ { 1 d ^HK» < 33 > 



d r P <9£ 

, ! / / ; - 

while we obtain for the Kondo order parameter A 



4 i-oo ~~JDee E 



4f3J k X{l - t^L / Dz— ±-= / Dee* 



1 r+PD l sinh VAL , , 

PDJ-pd cosh (*±*) + cosh (VA) V ^A V 7 

The numerical solution of the saddle point equations as a function of j 
and ^j provides us with the phase diagram in Fig.l, that we discuss in the 
next section. 



3. Conclusions 

We study in this paper the phase transitions in a system represented by a 
Hamiltonian that couples the localized spins of a Kondo lattice [||, |I"3||with 
random, long range interactions, like in the SK model for a spin glass |T0|. 



Using functional integrals techniques and a static, replica symmetric 
Ansatz for the Kondo and spin glass order parameters, we derive a mean 
field expression for the free energy and the saddle point equations for the 
order parameters. The Kondo and spin glass transitions are both described 
with the mean field like static ansatz that reproduces good results in two well 
known limits: when 4 = we recover the exact solution for the quantum 



Ising spin glass [14, 18| while for J = we recover the mean field results 



for the Kondo lattice |3], |13|.The use of the static ansatz is justified at finite 
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temperatures ||15||. Numerical solution of the saddle-point equations allow us 



to draw the magnetic phase diagram in the JkVsT plane, for fixed value of 
J, presented in Fig.l. 

Figure 1 shows three different phases. At high temperatures, the "nor- 
mal" phase is paramagnetic with vanishing Kondo and spin glass order pa- 
rameters, i.e. A = q = 0. When temperature is lowered, for not too large 
values of the ratio Jr/J, a second-order transition line is found at T = T$g 
to a spin glass phase with q > and A = 0. Finally, for large values of the 
ratio Jr/J, we recover the "Kondo" phase with a non-zero A value and q = 0: 
the transition line from the paramagnetic phase to the Kondo phase for tem- 
peratures larger than Tgc is a second-order one and occurs at a temperature 
very close to the one- impurity Kondo temperature Tr. On the other hand, 
the transition line from the spin-glass phase to the Kondo phase, for tem- 
peratures smaller than T$g, is a first-order one and it ends at J K at T = 0. 
When the temperature is lowered, the transition temperature does not vary 
very much with the value of Jr/J] the separation between the spin-glass and 
the Kondo phases departs completely from the behaviour of Tr and looks 
like the separation between the magnetic and Kondo phases when these two 



phases are considered |21||. We can also remark that we get here only "pure" 
Kondo or SG phases and never a mixed SG-Kondo phase with the two or- 
der parameters different from zero; this result is probably connected to the 
approximations used here to treat the starting Hamiltonian. 

The diagram shown in figure 1 can explain the magnetic phase diagram 
observed above the Curie temperature for the CeNi\- x Cu x || for small x 
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values when there is a transition from a spin-glass state to a Kondo state 
and then to the intermediate valence compound CeNi; however, there is no 
experimental information on the precise nature of the SG-Kondo transition 
and our model cannot be checked from that point of view. There is also 
probably a SG-Kondo transition in the CeCoGe?,_ x Si x alloys, but there the 
experimental situation is even more complicate than in the preceding case. 
An unsolved basic question concerns also the existence or not of a "mixed" 
SG-Kondo phase in Cerium disordered alloys and this problem is worth of 
being studied experimentally in more detail. Thus, further experimental work 
is necessary, but our model yields a new striking point in the behaviour of 
heavy fermion disordered alloys in the vicinity of the quantum critical point. 

Appendix 

We present here a detailed derivation of the main equations of the paper. By 
introducing eqs. (]l|)-(0) in e q- © we obtain for the s-d exchange part of the 
action: 

A K = -pj k j2 E Ev>l(w)Vw(w' + n)pL(wOM™-«) (35) 

i w,w',Cl a 

where w = (2n + l)ir and Q = 2mr. In the mean field spirit we want to 
introduce the spatially uniform and static Kondo order parameter in eq. 
fll4|), then we take Q = 0, re-order the operators and separate the sites in 
eq.(3D) with the introduction of a iV _1 factor, which gives: 



A * « +tEEwKWEW^M. (36) 
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that is eg. (|TT|) . 

We find it convenient to introduce the Kondo order parameters in eq.([l4|) 
by means of the identity 

/oo 
-oo 

x[A CT iV -Y.vlH^H]e" JkNl ^ Xl+ ^\ (37) 

and using the integral representation of the 5-funtion: 

5( x -x ) = — due tu{x ~ xo) (38) 

Z7T J-oo 

we may write the partition function by combining eq.(|7|) to eg. (|T3"|) : 

z = jy^uKjy/-^ 

xe - Z eff (39) 

where 

z eff = [d(^ijj) I z%V) e ^o+A SG 



xe" j ' w j ' w (40) 

The saddle point values of X a , A^., u a , u\ are otained by extremizing the 
exponent in eq.(^): 

f3J k \- a = iu a 
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4 = tt £<$»M«0> 



iV 



.;■«' 



^ = iE(A(*H ( 41 ) 

where (...) = / D(ip^)D{i^ (p)e Aeff (...) from eg. (|4*0|) . Introducing the saddle 
point values of eq. flHf ) into eq. (|39|) , we obtain: 

Z = e- 2N( " J ^ x Z eff . (42) 

From eq. ( 10 ) , Z e // is now the integral of a quadratic form in the ip^ , if 
variables, that can be integrated out to give: 

Z eff = Z° d Z SG (43) 

where Z° is the partition function of the free conducting electrons: 

ln(Z d °) = 2^1og|7^H|, (44) 

and 

lij 1 ^) = iw $ij ~ PUj, (45) 

is the inverse Green's function for the d-electrons. The quantity Zsg hi 
eq.fl4"3]) is the partition function for the localized f-electrons: 

Y,Y,9i j 1 Mijl(w)ij ja ( y w) + A SG 
Z SG = / D{^)e w id (46) 



where the inverse Green's function for the localized, non-interacting f-electrons 
is now modified by the Kondo interaction: 

gj(w) = (iw - Pe )6ij - (3 2 J 2 k ^\ ll3 (w) (47) 
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and Asg is given in eg. (|I~2j) . 

The interesting part of the free energy is given by: 

/^=-^<<ln(|j)}> (48) 

where the double bracket indicates a configurational average over the random 
variables Jij, with the distribution probability in eq.([llj). Using the replica 
method we obtain from eq. fl42|) and eq.( 



(3F = 2(3J k \^\ - hm^-[Z n (SG) - 1] (49) 

where a = l...n is the replica index and: 

/n 
UDUtM exp {£ £ 9^(w) £ t(ui)W«;)} 
a. wa lj a. 

xJI((e ° )) (50) 

The operators S ia are bilinear combinations of i/j' ia . a (w), ?/w(u>) from 
eg. (|T3|) , then after performing the average in eq.(|50|) we must use standard 
manipulations with gaussian identities [16,17] to linearize the exponent in 
eg. (|50|) . We obtain: 

Z n (SG)= Hdq aP e "J A({q a p}) (51) 

where 



» 3 <Ta\W) 



HUa/s}) = / Dto)exp{^ ^ ft/(w)^L(w)^ 

i j w,<j,a 

+P 2 J 2 Y,Q*0H S A} (52) 

a,/3 i 
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We obtain for Z n (SG) at the replica symmetric saddle point: 

q aa = q + X (53) 

Z n (SG) « e -^ J)2N ^ 2+2q ^A(q, x) (54) 

where 

/OO /"OO 

n Dz i n / n j%*j«(?> %> fe}> &* » (55) 
-00 . a J -co ■ 

and 
I a (q, X, {zj}, {taj}) = J D{$J) a ) exp {J2 J2 G ^LHVi a (w)Vw( w )} 

(56) 

where 

G--L(^) = 9v\w) - S ij a[ y /2q(3Jz j + ^(3J£ aj \. (57) 

Introducing Eq.(|56|) and Eq.(|55|) in to Eg. ([541), we obtain from Eq. (|49|) 
the expression for the free energy in Eq.fllTSD of the main text. 



w,cr 



> CT hj 



Y, ln \ G ijla{w) 
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Figure captions 



Figure l:Phase diagram in the T — Jk plane as a function of T/J and Jk/J 
for fixed J = 0.05/}, where D is the conduction bandwidth. The dotted line 
represents the "pure" Kondo temperature T^. 



21 







f/l 



